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Introduction. 



The purpose of this paper is to present an extended notion of strong solution to S.D.E.'s 
driven by Wiener processes. These solutions can be defined on rather general spaces, in the 
context of Dirichlet forms. 

More interestingly, they are not always given by flows of maps but by flows of Markovian 
kernels, which means splitting can occur. Coalescent flows also appear as solutions of these 
S.D.E.'s. Conditions are given under which coalescence and splitting occur or not. 

A variety of examples are studied. The case of isotropic Sobolev flows on the sphere or 
on the Euclidean space shows in particular that splitting is related to hyp erinst ability and 
coalescence to hyperstability. These notions (which will be developped in sections 9 and 10) 
are related to the explosion of the Lyapunov exponent toward +00 and —00. 

The typical example we have in mind is the Brownian motion on a Riemannian manifold. 
We consider a covariance on vector fields which induces the Riemannian metric on each 
tangent space. When the covariance function has enough regularity, it is known that one can 
solve the linear S.D.E. driven by the canonical Wiener process associated to this covariance 
(or in other terms to the local characteristics associated to this covariance (see section 3 
below)) and get a multiplicative Brownian motion on the diffeomorphism group, which moves 
every point as a Brownian motion (see Le Jan and Watanabe [^3| or Kunita [fLSfl). But models 
related to turbulence theory produce natural examples where the regularity condition is not 
satisfied. Except for the work of Darling |7]], where strong solutions are not considered, these 
S.D.E.'s have not been really studied. The idea is to define the solutions by their Wiener 
chaos expansion in terms of the heat semigroup. We call it the statistical solution. A similar 



expansion was given by Krylov and Veretennikov in |T7|], for S.D.E.'s with strong solutions. 

In this form, they appear as a semigroup of operators, and the fact that these operators 
are Markovian is not clearly visible in the formula. To prove this, we consider an independent 
realization of the Brownian motion on the manifold and couple it with the given Wiener pro- 
cess on vector fields using certain martingales. Then, the Markovian random operators which 
constitute the strong solution are obtained by filtering the Brownian motion with respect to 
this Wiener process. They determine the law of a canonical weak solution of the equation 
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given the Wiener process on vector fields. This construction has been adequately generalized 
to be presented in the case of symmetric diffusions on a locally compact metric space. It is 
a convenient and well studied framework but this assumption could clearly be relaxed (in 
particular to the framework of coercive forms). Relations with particle representations and 
filtering of S.P.D.E.'s can be observed (see Kurtz and Xiong |19|). 

The example of Sobolev flows is studied in details on Euclidean spaces and spheres and 
is of major interest especially in dimension 2 and 3 where an interesting phase diagram 
is given in terms of the two parameters determining the Sobolev norm on vector fields : 
The differentiability index and the relative weight of gradients and divergence free fields 
(compressibility) . 



Some of these results have been given in the note |2l| and a preliminary version of this 
work was released in |[22|| . They are directly connected and partially motivated by a series 
of works of Gawedzki, Kupiainen and al on turbulent advection (0, and |H}|) 



Acknowledgment : we wish to thank an anonymous referee for his careful reading and 
his suggestions which helped to improve the manuscript in many ways. 

1 Covariance function on a manifold. 

Let X be a manifold. A covariance function C on T*X is a map from T*X 2 in R such that, 
for any (x, y) G X 2 , C restricted to T*X x T*X is bilinear and such that for any n-uples 
(£i,...,£n) of T*X, 

For any £ = (x, u) G T*X, let be the vector field such that for any £' = (y, v) G T* X, 

<c € (v),*> = c(e,o- 

Let H be the vector space generated by the vector fields C^. Let us define the bilinear 
form on H , (.,.)# such that 

(Q,CV)h = C(£,0- (1.2) 

As equation (1.1) is satisfied, the bilinear form (.,.)# is a scalar product on H . We denote 
the norm associated to (., .)#. 



Let H be the separate completion of H with respect to ||.||ff. (H, (■, -)h) is a separable 
Hilbert space and we will design it as the self-reproducing space associated to the covariance 
function C . H is constituted of vector fields on X and for any h G H and any £ = (x, u) G 
T*X, 

{C i ,h) H = (h(x),u). (1.3) 

Let (e k ) k be an orthonormal basis of H, then equation (1.3) implies that for any £ = 
(x,u) G T*X, 

Cg = y^ j {e k (x),u)e k . (1.4) 

k 

This equation implies that for any £' = (w,t>) G T*X, 

it 

Therefore 

C = ^e fe (g)e fc . (1.6) 
fe 

Remark 1.1 On i/ie ot/ier hand, if we start with a countable family of vector fields (V^)*., 
sttc/i that for any £ = (x,m) G T*X ; (Vfc(a:), w) 2 < oo, it is possible to define a covariance 
function on X by the formula 

c = Y J v k ®v k . 

k 

Examples of isotropic covariances are given in section 9 and 10. See also Q. 
Assume now a Riemannian metric (., .) x is given on X, the linear bundles TX and T*X 
can be identified. We will now suppose that the covariance is bounded by the metric i.e that 

for any £ = (x,u) G T*X. Note that this condition implies that |/i(x)| x . < \\h\\H for any 
heH. 

Let us denote by m(dx) the volume element on X. Given any different iable function / 
such that |V/| is square integrable, we can map it linearly into Df in the Hilbert tensor 
product L 2 (m)® H setting (Df,g<g)h) = f x g(x)(Vf(x), h(x)) x m(dx) for all g G L 2 {m) and 

heH. 
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Note that 

\\Df\\ 2 H (x) < |V/(x)| 2 (1.7) 

(This notation comes from the identification L 2 {m)®H with the L 2 space of H valued func- 
tions on X) and that 

\\Df\\h {m)m < J \Vf\ 2 dm. (1.8) 

2 Covariance function bounded by a Dirichlet form. 

We can extend these notions to the framework of local Dirichlet forms. Let X be a lo- 
cally compact separable metric space and m be a positive Radon measure on X such that 
Supp[m] = X. 

Let (£, J 7 ) be a regular Dirichlet space, T C L 2 (X, m). We will suppose that the Dirichlet 
form is local and conservative. We will denote Pt, the associated Markovian semigroup, A its 
generator and T>(A) the domain of A. We will also suppose that m is an invariant measure, 
hence that P t l = 1. We will also assume that for any / G jF b = L°°(m) n J 7 , there exists 
r(/, /) G L x (m) such that for any g G 

2£(/0, /) - S(f, g) = J gT(f : f) dm. (2.1) 

T can be extended to T and we denote r(/, g) the L 1 (m)-valued bilinear form on J 72 , where 
for any (/, g) G J 72 , V(f, g) = \ (T(f + g, f + g) -F(f - g, f - g)). A sufficient condition 
for the existence of V (see corollary 4.2.3 in [Q) is that T>(A) contains a subspace E of T>(A), 
dense in J 7 , such that 

V/ EE, f 2 G V(A). 

Then, for (f,g)eE 2 , 

T(f,g)=A(fg)-f Ag-gAf. (2.2) 

A necessary and sufficient condition for the existence of the energy density (or carre du 
champ operator) V is given in theorem 4.2.2 in [j3|. 



5 



Fundamental example 2.1 X is a Riemannian manifold with the metric (.,.), m is the 
volume measure, T = i/ 1 (X) and for any (f,g) G T 2 , 

£(f,g) = \ [ (v/,v ff > dm. 

In this case, T(f,g) = (V/, V#). 

Let if be a separable Hilbert space and Z) a linear map from JF into the Hilbert tensor 
product L 2 {m)®H such that, for any / G T 

\\Df(x)\\ 2 H <T(f,f)(x) (2.3) 

m(dx)-a.e. The most interesting case is when there is equality in equation (2.3). 

We define the covariance function C as a bilinear map from T x T into L 2 (m <g> to) by 



{C(f,g),u® v) L 2 {m<S)m) = / (Df(x),Dg(y)) H u(x)v(y) m(dx)m(dy). (2.4) 

Note that 

(C(/,/), M ® M ) L2(m0m) <2£(/,/)Hi 2(m) . (2.5) 
We will say that C is a covariance function bounded by the Dirichlet form 

Remark 2.2 Alternatively, we could define the covariance as a positive bilinear map from 
T ' x T in L 2 {m <g> to) (i.e such that for any Ui G L 2 {m) and any fi G T , 

[J2 U *® u i C(/<> /i) <*m® 2 > 0-) ( 2 - 6 ) 

such that 

(C(f,f),u®u) LHm ® m) < 2S(f,f)\\u\\ 2 L2{m) (2.7) 

and construct as before a Hilbert space such that (2.3) holds. 

Indeed, C induces a linear map C from L 2 {m) £g> T into L 2 {m)®T such that 

{C(u® f),v®g) L 2 {m) ® T = (C(f,g),u®v) L 2 {m ® m) . 

We define H as the separable closure of the space H spanned by elements of the form 
C(u<8) f), with u <S> f G L 2 {m) (g) T , and equipped with the scalar product 

{C{u® f),C(v <g> g)) H = (C{f,g),u®v) L 2 {m ® m) . 
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And for f G T , Df is defined such that for any u <g> v <S> g G L 2 (m) <g> L 2 {m) <g> T , 
(Df,u®C(v®g) > L 2 (m) ® H = (C(f,g),u®v) L 2 {m ® m) . 

For any h G H and / G define D^/ = (-D/, /i)h which belongs to L 2 (m). Then for any 
orthonormal basis (ek)k of if, 

C7 = X)^®I>e fc . (2.8) 

Moreover, for any / G J 7 , 

\\Df\\ 2 H = J2( D eJ) 2 - (2-9) 

k 

Remark that condition (2.3) implies that for any finite family (u,, /,) G L°°(m) x JF, 

J^UiUjDifiJj) < J^UiUjFifiJj), (2.10) 

where D(f,g) denotes {Df(x),Dg(x))u = Ylk-De k fi x )D ek g(x). When the tij are step func- 
tions with discontinuities in a set of zero measure, (2.10) is satisfied as £V . u i UjD(f i , fj) = 
\D(J2iUifi)\ 2 . Then, we can extend to any family (tij) by density in L 2 (T(f, f)dm) for every 

Remark 2.3 f£ is c/ear t/iat <jwen a covariance C on T*X as in Section 1, we can build the 
self -reproducing space H consisting of vector fields and the mapping D : ii" 1 (X) — > L 2 {m)®H 
so as to construct a covariance function as in Section 2. Now suppose conversely that we 
have a separable Hilbert space H , a linear map D and a covariance C as in Section 2, and 
suppose we are in the Riemannian case. The condition \\Df(x)\\ 2 H < F(f,f)(x) = |V/(a;)| 2 
implies that C(f,g)(x,y) depends only on V/(x) and Vg{y), and so there is a covariance 
C say on T*X so that C(f,g)(x,y) = C(Vf(x),Vg(y)). So in the Riemannian case, any 
Section 2 covariance function reduces to a Section 1 covariance function. 

Further, we can now assume without any loss of generality that the separable Hilbert space 
H is the self -reproducing space corresponding to C and thus consists of vector fields. 

Remark 2.4 The bilinear mapping D is a derivation : for any h G H and any f G T such 
that f 2 G T, 

D h f 2 = 2fD h f. (2.11) 
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Note that in the Riemanian manifold case (fundamental example 2.1), Dhf = V^/ when 

r = d. 



Proof. We first make the remark that 

J2( D e k f - VDejf = D(f, f) - 4fD(f, f) + 4/ 2 D(/, /). 

k 

Integrating this relation with respect to m and using (2.10), we get that 

[j2(D e J 2 -2fD e J) 2 dm< f (T(f,f)-AfT(f,f)+AfT(fJ)) dm = 0. 

This implies that for every k, D £k f 2 — 2fD £k f = 0. □ 

3 Construction of the statistical solutions. 

In the fundamental example 2.1, when X is a Riemannian manifold, C is smooth and when 
equality holds in (2.3), it is well known (see |23 and that a stochastic flow of diffeomor- 



phisms on X can be associated with C. Then, with the notations of definition 2.1 in |[23|| . 
the local characteristics of the flow are (A, L) where A = C and L is the Laplacian on X. 

The object of this section is to show that in the general situation considered above, it is 
always possible to define a flow of Markovian kernels associated with C and {£,3-") (which 
is induced by the stochastic flow when C is smooth). 

Let be given a covariance function C bounded by a Dirichlet form (£,J-) on a locally 
compact separable metric space as in the preceding section (equation (2.3) is satisfied). Let 
W t be a cylindrical Brownian motion on H defined on some probability space (Q,A,P), 
i.e a Gaussian process indexed by H x R + with covariance matrix cov(Wt(h),W s (h')) = 
sAt (h, h')n- Set = W t (ek). (W^; k G N) is a sequence of independent Wiener processes 
and we can represent W t by Y2k ^t e k- Informally, the law of W t is given by 

z 

Let T t = o-{W*; k e N; s < t) = a(W s ; s <t). 



Proposition 3.1 Let S® = P t . We can define a sequence S™ of random operators on L 2 (m) 
such that E[(Stf) 2 ] < P t f 2 in L l {m) and 5" is T t -measurable, by the recurrence formula, 
in L 2 {m ® P) (i.e in the Hilbert tensor product L 2 (m)(£)L 2 (P) ) 

S? +1 f = Ptf + J2 f S s( D e k Pt-sf)) dWl (3.1) 

k J ° 

Remark. The stochasic integral in equation (3.1) here makes sense as a Hilbert valued Ito 
integral. Recall that given a real Wiener process W t and a Hilbert space H, for any F progres- 
sively measurable in L 2 (P W <g) dt)®H and any h E H, (J* F(u)dW u , h)n = fs(P( u )i h)ndW u 
and E[\\ f s F(u)dW4%} = fj\F(u)\\ 2 H du. 

Proof. Suppose we are given S", a .^-measurable random operator on L 2 (m) such that 
E[(S?f) 2 ] < P t f 2 . 

Let / G L 2 {m). For any positive t, P t f G T and D ek P t ^ s f is well defined and belongs to 
L 2 (m). 

E[(Sr l ff] = Wf + Y. f ' E [(S?(D ek P t _J)) 2 } ds m-a.e. 

< (P t f) 2 + f P s {\DP t _ s f\ 2 ) ds 

Jo 

< (P t f) 2 + [ P s (r(P t - s f,p t - s f))ds. 

Jo 

For / G L°°(m) fl L 2 (m), £ P 3 ((P t _ J) 2 ) = P s (r(P t _J, P t _ s f)) and 

Ptf = (Ptf) 2 + f P s (r(P t - s f, P t - S f)) ds. (3.2) 
Jo 

An approximation by truncation shows that equation (3.2) remains true for / G L 2 (m) and 
E[(Sl l+1 f) 2 ] < PJ 2 . □ 

Remark. The definition of S™ is independent of the choice of the basis on H. 

In the following, we will use the canonical realization of the processes W^. They will be 
defined as coordinate functions on Q = C(IR + ,IR) N , with the product Wiener measure P. 
We note 9 t the natural shift on Q, such that W^ +s - = o 9 t . 

Recall that an operator on L 2 (m) is called Markovian if and only if it preserves positivity 
and maps 1 into 1 (or more precisely, if m is not finite, if its natural extension to positive 
functions maps 1 into 1). 
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Theorem 3.2 The family of random operators 5" converges in L 2 (P) towards a one pa- 
rameter family of T t -adapted Markovian operators S t such that 

a) S t+S = S t (S s o 6 t ), for any s,t > 0; 

b) V/ E L 2 {m), S t f is uniformly continuous with respect to t in L 2 (m <g> P); 

c) E[(S t f) 2 ] < P t f 2 , for any f E L 2 (m); 

d) S t f = P t f + £ fe JJ S s {D ek Pt- s f) dW s k , for any f E L 2 (m); 

e) S t f = f + £ fe £ S s (D e J) dW k s + JJ S s (Af) ds, for any f E V(A). 

S t is uniquely characterized by c) and d) or by a), c) and e). When Y = D, we call it the 
statistical solution of the S.D.E. (see 3.22 below) 

V/ E V{A) : df(X t ) = Y,D e J{X t ) dW? + Af(X t ) dt. (3.3) 

k 

Note that this S.D.E. does not always have a strong solution in the usual sense. 

Proof. The convergence of 5" is immediate since for any n > 1, J"/ = S™ f — S™~ 1 f is in the 
Hilbert tensor product of the n-th Wiener chaos of L 2 (P) with L 2 (m), Stf = Ptf + Y^=i J-tf 
and (P t f) 2 + £n>i£pf/) 2 ] = linw,£[(S t "/) 2 ] < P t f 2 . It is clear that S t is ^-adapted 
and satisfies c). d) is obtained taking the limit in the recurrence formula of the proposition. 
Since 

Jtf= E / P sl D ek P S2 _ sl ...D ekn P t - s J dW£...dW*: 

U. U J0<S 1 <...<Sn<t 

we have J" +s = £ fc<n J k { J™~ k ° t ) (the k-th term corresponds to 

E / P sl D ek P S2 - sl . . . D ekn P t+s _ s J dW^...dW k : ). 

ki,...,k n ^ 0<s 1 <...<s k <s<s k+1 <...<s„<t+s 

We deduce a) from this relation. 

The uniqueness of a solution of d) verifying c) follows directly from the uniqueness of the 
Wiener chaos decomposition, obtained by iteration of d) : Let T t design another solution of 
d) and c) then for any / E L 2 (m) and any integer n, 

T t f = Sr l f+ E / T Sl D eki P S2 _ Sl ...D ekn P t - s JdW* l \..dW* : . 

k, k J0<S!<. ..<Sn<t 
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The second term of the right hand side of the preceding equation is orthogonal to the first 
one since its integrands are L 2 . Indeed : 



E E I (T Sl D eki P S2 - Sl ...D ekn P t _ s j) 2 d Sl ...ds n 

u U0<s 1 <...<s n <t v 7 

r»l j ■ ■ ■ ^rffh 

< I Ps 1 (\De kl P S2 -s 1 .-.D ekn P t „ s J\ 2 )ds 1 ...ds r 

k , z, >>0<si<...<s„<t 

< I P S2 (\D ek2 P S3 . S2 . . . D ekn P t . s J\ 2 ) ds 2 ...ds r 

fc 2 ,...,ifc n • / °< S 2< -< S ™<* 

using equation (2.3) and (3.2) and by induction is smaller than Ptf 2 . 

This proves that the Wiener chaos decomposition of T t f and S t f are the same and 
therefore T t = S t . 

Proof of b). Let us remark that for any positive e, S t+£ — S t — S t (S £ o t — I). As S t and 
S £ o 6 t are independent and m is invariant under P t , for any / G L 2 {m) 

j E[(S t+£ f - S t f) 2 ] dm < J E[P t (S £ o 6 t f - ff] dm 

< J E[(S £ o6J-f) 2 ] dm 

< J(P £ f 2 -2fPJ + f 2 )dm 

< 2||/|| L2(m) ||/-P £ /|| L2(m) . (3.4) 

Therefore, lim e ^ \\S t+£ f - S t f\\ L *(m®p) = °, uniformly in t. □ 

Remark 3.3 Note also the convergence in L 2 (m £g> P) of P £ Stf toward Stf when e — > 0. 
Indeed \\P £ S t f - S t f\\ 2 L2(m ^ p) = E[\\P £ S t f - S t /||| 2(m) ] and \\P £ SJ - S t f\\\ 2{m) converges 
towards when e goes to and is dominated by 4||S , t /||| 2 ( m )- 

Proof of e). Let us remark that for any e and t positive, 

St + ef-S t f = S t ^P £ f + J2fs u o9 t (D ek P £ _ u f)dW^o9 t -f^j 

= S t (PJ -/) + £ [ t+£ S s (D ek P t+£ _ s f) dWt (3.5) 
k Jt 
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Hence using (3.5) for t = H and e = H, for / G V(A), 

Stf-f-J2 f S s( D eJ) dW k s - f S s (Af)ds = 
i, Jo Jo 



n-l 

E 

i=0 



S± t (P*f - f) + E / S s (D ek Pi+i t _J) dW k s 



S.(Af) ds-J2 i S s (D e J) dW k s 



= Ai (n) + A 2 (n) + A 3 (n) , with 



n-l 



A 1 (n) 
A 2 (n) 
A 3 (n) 



J2S ±t (P±f-f--Af); 

f * n n fl 

1=0 

n-l ,!+l t 

/ " (S lt (Af)-S s (Af))ds; 

1=0 n l 

n-l .i±l 4 

EE A 5 s (D efc (p^ t _ s / - /)) dwt 



(3.6) 
(3.7) 
(3.8) 



=0 k 



First, using the fact that m is P t -invariant, 



\\Ai(n)\\ L 2 {m ® P) < n\\P±f - f - -Af\\ L 2 {m) = o(l). 



(3.9) 



After, we remark that 



i+lj 



{S u (Af)-S a (Af)) ds 



2 














n J± 




-M n 



i+lj 



\\S ±J (Af) ~ S s (Af)\\l Hm0P) ds. 



As St(Af) is uniformly continuous in L 2 {m®P), there exists such that lim^o^O^) = 
and \\S ±J (Af) ~ Ss(Af)\\l 2[mrm < e (i) for any s G [if, i±if]. Hence we get 



(S u (Af) - S s (Af)) ds 



L 2 {m®P) 



^ — £ - 

n 2 n 



and \\A 2 (n)\\ L 2 {m ® P) = o(l). 

At last, as the different terms in the sum in equation (3.8) are orthogonal, 



n-l 



\\A 3 (n^ 2 



L 2 (m(g.P) 



= EE * 



i=0 k 



S s {D ek {P^ t _J - f))dW t 



dm 
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n-l 



<- E 



i=0 




< n 




\D(P^ t _J - f)\ 2 ds dm 



\D(P s f - f)\ 2 dm ds 



< n / £(P s f-f,P s f-f) ds. 
Jo 

As \im s ^ S(P s f - f,PJ - f) = 0, ||A 3 (n)|| L 2 (m8P) = o(l). 

Taking the limit as n goes to oo, this shows that \\S t f — f — J2k lo S s {D ek f) dW k — 
f*S s (Af) ds\\ L 2 {m ® P} = 0. □ 

Proof that a), c) and e) imply d). Take / G L 2 {m) and e positive, assuming e), 

S t PJ - P t PJ -Vf S s (D ek P t _ s PJ) dW k s = 
k Jo 

54ti t (P t _iti t Pe/) - 5 lt (P t _ ±t P e /) - J] / " S s (D ek (P t _ s P £ f)) dW k s 



n-l 

E 

7=0 



Bi (n) + 5 2 (n) + 5 3 (n) , with 
fii(n) 



71-1 

EE A " 5 s p efc (P t _^ t P £ /-P t _ s P £ /))^ s fe ; 
i=o fc J 



(3.10) 



71-1 



i=0 



S 2 (n) = -J>± t [Pt-HPef ~ Pt-^tPef ~ -AP^PJ ; (3.11) 



B z{n) = ^ (S s - S± t )(AP t _i+i t P £ f) ds, since 

^— » / i . n n 

„_n •/ —I 



(3.12) 



Si±i t (P t _i±i t PJ) 



i=0 

5i t (P t _i±i t P £ /) 



n 



^(AP^^PJ) <fe 



+E 



SsD^P^PJ dW k s 



Since the different terms in the sum in equation (3.10) are orthogonal, 



\BAn) 



\L 2 {m®P) — 



n-l 

E 

7=0 




|P>(P t _l±l t P £ / - P- S P £ /)| 2 <fe dm 



it 



n-l /.i+ly 



< V/ £(P t _±ti t P e f -P t . a P e f,P t _^ t P e f -P t - s PJ)ds 

7=0 n 1 

< n I" £(P s P £ f-P £ f,P s P £ f-P £ f)ds (3.13) 

./o 
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as £(Ptf, Ptf) < £(/, /) for any positive t and any / G L 2 (m). 

Equation (3.13) implies that \\B 1 (n)\\ L 2 {m ® P) = o(l) (as \im s ^ £(P s PJ - P £ f,P s P £ f - 
Pef)=0). 

\\B 2 (n)\\ L 2 {m ® P) < Y] \\S± t P t -±t P ef ~ P t -i±±t P ef ~ -AP t _i±± t P £ f \W{m®P) 

f * n \ n n 71 n/ v/ 

< n\\P ± PJ - PJ - -APJ\\ LHm) . 

n n 

Hence, ||S 2 (n) \\ L 2( m ®p) = o(l). 

Note that if Q t / = e) implies that for any / G V(A), 

Qtf = f+f Qs(Af)ds. 

Jo 

Then -jj-Q 8 P t _ s f = for any / G L 2 (m) and < s < t (then P t - S f G ©(A)) and we have 
Q t f = P t f. With this remark and the fact that a) and c) are satisfied, we see that b) and 
equation (3.4) are satisfied (see the proof of b). Using (3.4), we have 

\\(S s -S it )(AP t _^ t P £ f)\\l 2{m ^ P) < 

< 2\\AP t _ 1±lt PJ\\ LHm) \\AP t _ 1±lt PJ - P^AP^^PJU^) 

n n n n 

< 2\\AP £ f\\ LHm) \\AP £ f - P s ^ t AP £ f\\ LHm) 

< A\\AP £ f\\ 2 L2(rn) . 

Hence, 

n— 1 i+i j 

\\Bz{n)\\l^ P) < E^// \\( S s- S u)( AP t-^tPef)\\Urn®P)ds 

4/2 

< —11^/11^). 

Taking the limit as n goes to oo, this shows that d) is satisfied for P £ f, with / G L 2 {m) and 
e positive. 

At last, since \\S t P £ f - S t f\\ L 2 (m ® P) < \\P £ f - /|| L a (m) (because c) is satisfied), \\P t+£ f - 

Ptf\\mm®P) < \\Pef - f\\mm) and || £ fe £ S s (D ek P t _ s (P £ f - f))dWX Hm ®P) < t£(P e f - 
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/, P e f — f). Taking the limit when e goes to 0, we prove that d) is satisfied for any / G L 2 {m). 
□ 



Proof that St is Markovian. 

A more concise proof of this fact has been given in ||21|| , relying on Wiener exponentials 



and Girsanov formula. The advantage of the following proof is to be more explanatory, to 
give a relation with weak solutions and to yield a construction of the process law associated 
with the statistical solution St- 

Let (ft 1 , Q, Q t , X t , P x ) be a Hunt process associated to (£,-7 r ) (see [[□]]), we will take a 
canonical version with Q' = C(M + , X). Let Wt = W^ek be a cylindrical Brownian motion 
on H, independent of the Markov process X t . 

Let M. be the space of the martingales additive functionals, ^-adapted such that if 
M G M, E X [M?] < oo, E x [M t ] = q.e. and e(M) < oo where e(M) = sup t>0 j- t E m [M 2 } 



(with P m = J P x dm(x)). (Ai,e) is a Hilbert space (see |TT| ). 

For / G J 7 , M* G M. denotes the martingale part of the semi-martingale f{X t ) — f(Xo). 
For g G C K (X) C L 2 (T(f, /)tZm)fl, we note g.M f G M the martingale ftg(X a ) dM f s , then 
M = {Er=i^- M/! ! n e N, g t G C K (X), f\ G T} is dense in M (see lemma 5.6.3 in |TT[), 
and e(^2 i gi.M^ i ) = jj^ijl 9i9j^(fi, fj) dm (see theorem 5.2.3 and 5.6.1 in fLTf ). 



Lemma 3.4 For every (M,N) G M x 7W, t/iere exists r(M, AT) G L x (m) s^c/i t/iat 

(M,N) t = [ T(M,N)(X S ) ds } (3.14) 
where (., .) t zs t/ie usim/ martingale bracket. And for (f,g) G JF, r(M^,M 9 ) = T(f,g). 
Note that lemma 3.4 implies that e(M, iV) = ± / T(M, JV) dm. 

In the fundamental example 2.1, X t is the Brownian motion on X, M[ is the ltd integral 
J*(df(X s ), dX s ), T is the inverse Riemannian metric and M. can be identified with the space 
of 1-forms equipped with the L 2 -norm associated with the metric. 



Proof. When / G JF, it follows from theorem 5.2.3 in [] i] that 

-t 



(Mf,Mf) t = ! T(fJ)(X s ) ds. 
Jo 



1 Ck(X) design the space of functions continuous with compact support. 
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For M = £\ hi.M fi , N = £\ kj.M®, two martingales of M , 

(m, N) t = V / //,^rf/,/ // )(.vj d s = / r(M,iv)(x s ) d s , (3.15) 

with F(M,N) = ^2ijhikjT(fi,gj). T is a bilinear mapping from A4 x .Mo in L^ijn). Y is 
continuous since for any (M, AT) G A4o x .Mo, 

y |r(M,AT)| dm < J T(M,M)^T(N,N)^ dm 
< 2e(M)h(N)^. 

It follows that T can be extended to M. x Ai. 

Take MgM and an approximating sequence M n G -M . Then e(M n — M) converges 
towards 0, M n converges towards M in L 2 (P X ) and {M n ,M n ) t converges in L 1 (P ;r ) to- 
wards {M,M) t for almost every x (see section 5-2 in ||11|| ). This proves that (M,M) t = 
f*T(M,M)(X s ) ds. □ 

Lemma 3.5 If m is bounded, for any h G H , there exists a unique continuous martingale 
in M, N h such that for any f G T, e(N h ,Mf) = § f D h f dm and £ t (N h ,M f ) t = D h f(X t ). 
In addition, e(N h ) < \m(X)\\h\\ 2 and (N h ) t < \\h\\ 2 t. 

In the Riemannian manifold case (example 2.1), N t h = J*(h(X s ),dX s ) when T = D. 

Proof. For h = J2 k A^e^ G H, let us define a linear form, at on A4 such that for any 
M = 9i.M* G Mo, a h (M) = \ £? =1 / 9iD h fi dm. 

< ^\\h\\ 2 m(X)J2 j 9i9jD(fi,fi) dm< l -\\h\\ 2 m(X)e(M). 

i,3 

This proves that oth is continuous on .Mo and can be extended to a continuous linear form 
on M. such that oth(M) < ^\\h\\ a/ m(X)e(M). To this form is associated a unique N h G M. 
such that a h {M) = e(N h ,M). 

Note that for any g G CV(^0 and / G JF, we have J gDhf dm = f T(N h , g.M?) dm = 
J gT(N h , M*) dm. This is satisfied for every (7 G Ci<-(Jf), therefore for any / G JF, 

r(iv ft ,M/) = aj. 
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Note that we also have, for M G Mo 



r(N h , m) < \\h\\ r(M,M)5, 



which implies that (N h ) t < \\h\\H. □ 

Remark 3.6 When m is not bounded, N h can be defined as a local martingale such that for 
any compact K and any f G T , 1k-N h G M., e(lK-N h , M?) = \ f R Dhf dm. In addition, 
e(l K .N h ) < \m(K)\\h\\ 2 . 

Let 7h be a function on X such that J^(iV ei , N Ck ) t = 'Yki(Xt)- Lemma 3.5 implies that 
the matrix A = ((Ski — 7m)) is positive (as -^(N h ) t < |H| 2 )- Therefore, it is possible to find 
a matrix R such that R 2 = A. 

Remark 3.7 If for any f G T, V(f, /) = \\Dff H , then for any f G T h , 



D ek f = YliDeJ Iki(Xt) and the positive symmetric matrix P = ((jki)) ^ s a projector. In 
this case, R = I — P . 



Proof. Set Q{ = £ fc £ D eJ( X s) dN^, Q f G M, then for any M = £? =1 ft-M* G M , 




(3.16) 




This implies that for any M G M, e(Q f , M) = e(M f , M) and Q f = M f . 
Since by lemma 3.5, f t (M f , N ek ) t = D e J(X t ), we get that 

D e J(X t ) = ^ f (Q f ,N^) t = J2 D eJ(Xt)^-(N ei ,N ek ) t = [J2 D eJ 7« ) (**)■ 
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This relation implies that N^ k = J2i fo lki(X s ) dNp (this is easy to check, considering 
j- t (N e \M) t with M G Mo). From this, we see that lkl = Y.ilkilu (i.e P 2 = P).U 

Set W k = N e t k + Jo Rki{X s )dW l s and W t = £ fe W k e k . 

In the Riemannian manifold case, when r(/, /) = ||-D/||# for any / G denoting by 
C( Xjtt ) when u G T X X and £ = (x, w) we have : 

dW t = dW t + C(x u dx t ) - C(x t ,dWt(x t )) 

and dW k = dW t k + (e k (X t ),dX t )-J2(ek(X t ), ei (X t ))dWl. 

i 

In this case, R is a projector (see remark above). 
Lemma 3.8 (W k ) k is a sequence of independent Brownian motion. 
Proof. Since W k is a continuous martingale, we just have to compute ^(W k , Wj) t : 

±{W?,Wl) t = lkl + R 2 kl = 5 kl . 

This implies the lemma. □ 

Let \i be an initial distribution of the form hm, with h a positive function in L 2 (m)nL 1 (m) 
and for / G L 2 {m) define S t f by the conditional expectation 

S t f(X ) = E^[f(X t )\a(X ,W k - k G N; s < t)}. (3.17) 

(One checks easily that this definition does not depend on h.) Remark that as X t is Marko- 
vian and W t has independent increments, 

S t f(X ) = E^fiX^aiXo,^; k G N; s > 0)]. (3.18) 

In the same way, we see that S t satisfies the multiplicative cocycle property a). 

Lemma 3.9 For any f G T>(A) and \i an initial distribution absolutely continuous with 
respect to m, 

Stf = f + Y,f S s (D e J) dW k + f S s (Af) ds, P, a.s. 
k ^° ^° 
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Proof. For any / G T>(A), we have 

f(X t ) = f(X ) + M{ + f Af(X s ) ds. 

Jo 



(3.19) 



It is clear that E.[f*Af(X a ) ds\a(W k ; k G N; s < t)] = f* S s Af(X s ) ds, as (3.17) is 
satisfied. Let Z t = T, k Jo H s dW k G L 2 (a(W k ; keW, s<t)), 



E.[Z t M{\ = J2 E \f H k s d(W\M f ) s 
k Lio 

= J2 E \f tH s D eJ( X s)ds 
u Uo 



H k S s {D e J) ds 



= E 



Z t J2 f S s {D e J) dW\ 

u JO 



This proves that E\M t f \a(W k ; fceN; s<t)] = J2 k Jo S s D eJ dW k . □ 

Now, using uniqueness in theorem 3.2 and the isomorphism j between L 2 (cr(W k ; t > 
0; k G N)) and L 2 (a(W / t fc ; i > 0; fc G N)), we see that jS t = S u which implies that S t is 
Markovian. □ 



Proposition 3.10 For any f G T\>, the martingale 

p/ = Mf-J2 f D e J{X s )dW k 

v. Jo 



(3.20) 



is orthogonal to the family of martingales {W k ; k G N} 7 in the sense of the martingale 
bracket (i.e for any k, (Pf,W k ). =0). And for any (f,g) G T 2 , 

rt 



(P f ,P 9 )t = [ (r(f,g)(X s )-D(f,g)(X s )) ds. 
Jo 



(3.21) 



Proof. We just have to show that (pf, W k ) t = for every / G Tb and every k G N which 
is true as 

{Mf,W k ) t = (Mf,N e *) t = ! D e J{X s ) ds. 

Jo 



19 



Let (f,g) G J^ 2 , then 

(P f ,P 9 )t = (P f ,M°) t 

= (Mf,M°) t -J2[ D e J(X s )d(W k ,M°) s 
k Jo 

r(f,g)(X s ) ds-T [ D e J(X s )D ek g(X s ) ds 
k Jo 

r(f,g)(X s ) ds- [ D(f,g)(X s ) ds. U 
Jo 

Remark 3.11 In the case T(f,f) = \\Df\\ 2 H for any f G T, proposition 3.10 implies that 
P/ = and that 

M / = E f DeJ{X s )dWl 

k J ° 

From this, we see that the diffusion X t satisfies the S.D.E. 

f(X t ) - f(X ) = f D e J(X s ) dW k s + f Af(X s ) ds (3.22) 
t Jo Jo 

for every f G T>(A). Therefore (X t , W t ) appears as a weak solution of this S.D.E. and S t is 
defined by filtering X t with respect to W t . 

Let P x fij{duj') be the conditional law of the diffusion X t , given Xo and {Wt; t G M + } (it 
is independent of the choice of the initial distribution). Using the identity in law between 
W and W, we get a family of conditional probabilities P x ^(duu') on C(1R + ,X) defined m®P 
a.e. 

Remark that (with X t (d) = u'(t)) 

S t f(x,u) = J f(X t (u')) P x ,u{dw') m®P a.s. (3.23) 

Under P XjLU (duj')P(duj), X t (uo') verifies the S.D.E. (3.3). It is a canonical weak solution of the 
S.D.E. (3.3) on a canonical extension of the probability space on which W is defined. St is 
obtained by filtering X t with respect to W. 
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4 The n-point motion. 

Let P^ be the family of operators on L°°(m® n ) such that, for any (fi)i<i< n G L°°(m), 

Pi n) fi® ®f n = E[Sth® - ®$/ n ]. (4.1) 

p/ n ^ is a Markovian semigroup on L°°(m® n ) as St is Markovian and satisfies a) in theorem 
3.2. It is easy to check that P^ maps tensor products of L 2 (m) functions into L 2 (m m ). 

Proposition 4.1 For any family of probability laws on X absolutely continuous with respect 
to m, (fXi, 1 < i < n), 

Pfcl,»M^ <K) = [ P(du) ®u P,UM) (4-2) 

Jn 

defines a Markov process on X n (with initial distribution Cg)™^,) associated with P^ n \ We 
shall call this Markov process on X n the n-point motion. 

Proof. For every family of functions in L°°(m), (/$) i<i<„, m® n ®P a.e (with X^u^) = tu-(t)) 

n 

Sf n h® ...<g> f n (x 1 ,...,x n ,uj) = Y[S t fi(xi,u) 

i=i 

/n 
\[f\(XM))&UP^m)- (4.3) 
i=l 

We get the result by integrating both members of (4.3) with respect to P{duj). □ 

Let .D (n) be the linear map from H x T® n in L 2 (m® n ) such that for any (fi)i<i< n G T 
and h E H, 

n 

D^fi ® ... ® /„ = ^ /i ® ... ® ® ... ® /n- (4-4) 

i=l 

Proposition 4.2 For any (/;)i<;<„ G fl L°°(m), 

Sf n h ® ... ® / B = A ® ... ® / n + X T Sf n {D^h ® ... ® /„) dW? 

it 

./O 
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where 

n 

A^h®...®f n = Y,h®-® A h®-®fn 

i=l 

+ Yl X^ 1 ® -® D eJi® ...®D e J j ® ...<g> f n . 

l<i<j<n k 

Remark. 1) For n — 2, the formula extends to functions in V(A) and A^ f <g> g — Af ® 
g + f®Ag + C(f, g), where (/, g) E (V(A)) 2 . 

2) Taking the expectation, we see that A^ is the infintesimal generator of P^ on 

(V(A) nL°°H) 82 . 

3) The formula extends to C 2 K (X n ) in the Riemannian manifold case (using for example 
the uniform density of sums of product functions and the regularizing effect of P® n ). 

Proof. This is just a straightforward application of Ito's formula applied to S t fi <S> ■■■ <S> S t f n , 
using the differential form of the equation satisfied by S t , e) in theorem 3.2. Taking the 
expectation and differentiating with respect to t, we get 

4 pi n) h®-®fn = i E\srh®...®f n ] 

at\t=o at\t=o 

= A (n) /i ® ... <g>/„. □ 



Remark 4.3 In general, m® n is not invariant under . 

5 Measure preserving case. 

We say that the statistical solution S t is measure preserving if and only if mS t = m a.s for 
all t (i.e m is invariant for St). When m(X) = oo, we use the natural extension of St to 
L 1 (m) or to positive functions defined m-a.e. 

Let us denote by Tk the set of functions of T which have compact support. 

Proposition 5.1 St is measure preserving if and only if J C(f,g) dm® 2 vanishes for all f, 
g in Tk- Moreover, define r t on L 2 {TA by Wj? o r t = W^_ s — . Then the adjoint of S t in 
L 2 {m) is S* t =S t or t . 
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Remark 5.2 a) When f G T K , C(fJ) G L^m 02 ). 

b) In the Riemannian manifold case, the condition that J C(f,g) dm® 2 vanishes for all 
f , g in Tk is equivalent to assume that W t is divergent free in the weak sense, i.e that for 
any f G T K , j(W t ,Vf) dm = 0. (It follows from the identity E (J{W t , V/) dm) 2 = 
tfC(f,f) dm**.) 

Lemma 5.3 Assume that J C(f, g) dm® 2 vanishes for all f , g in Tk, then for every h G H , 
f, g in T, 

gD h f dm = - I fD h g dm. (5.1) 



Proof. For every h G H, (g, f)*—*f gDhf dm is a continuous bilinear form on T x T since 

ll^/lli 2M <^(/,/)II^IIL M - 

Take /, g in Tk H L°°(m) then f g G Tk (as the bounded functions of a Dirichlet space 
form an algebra) and, since D ek is a derivation, D ek (fg) = gD ek f + fD ek g. Using this 
property, we get 

J2 (/ (9DeJ + fD ek g) dm)) = E (/ D ^fs) dm ) 

= Jc(fg,fg)dm® 2 = 0. 

This implies that for every k, J gD 6k f dm = — J fD Ek g dm. To conclude we observe that 
both members of (5.1) are continuous in / and g and that Tk H L°°(m) is dense in T (since 



the Dirichlet form is regular, see section 1.1 in [p| ). □ 

Proof of proposition 5.1. Assume J C(f,g) dm® 2 = holds for every / and g in Tk- 
Let us remark that the expression of the n-th chaos of S t f is given by the expression 

Jtf= [ E P Sl D eki P S2 - Sl D ek2 ...D ekn P t _ s J dW^...dW^. (5.2) 

J0<s 1 <S2<-<s n <t ^ h 

From this expression, using lemma 5.3 and the fact that Pt is self-adjoint in L 2 (m), we get 
that for / and g in L 2 (m), 

j gJ r t l f dm = (5.3) 
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Making the change of variable u n -i + i = t — s i; we get that the adjoint of J™ is given by 
(J n Tt9= [ E Pu 1 D eki P U2 - Ul D ek2 ...D ekn P^ Un gdW^or t ..JW^or^ (5.4) 

J0<Ui<U'2<...<U n <t h 7, 

From this it is easy to see that S^g = o r t )g (as they have the same chaos expansion). 

Notice that S^l = 1. A priori the constant functions are not in L 2 (m), but there exists 
an increasing sequence in L 2 (m), g n such that g n converges towards 1. For any nonnegative 
function / G L 2 (m), 

J S t fg n dm = J fS* t g n dm. (5.5) 
This equation implies, taking the limit as n goes to oo, that 

mS t (f) = J fS;idm = m(f). (5.6) 

And we get that mS t — m a.s. Which ends the first part of the proof. 
Conversely, it follows from proposition 4.2 that for all /, g in V(A), 

S? 2 f®g-S t f®g-f®S t g + f®g- f Sf 2 C(f,g) ds 

Jo 

is a square integrable martingale. This result extends to /, g in T. Taking /, g in Tk, 
integrating with respect to m® 2 and taking expectation, we get that J C(f, g) dm® 2 vanishes. 
□ 

Remark 5.4 When S t is measure preserving, Pj^ is self-adjoint in L 2 (m® n ) and in par- 
ticular m® n is invariant under P^ n \ The associated local Dirichlet form is such that 

£ {2) (f ®9,f®9) = S(f, /) |b||i 2(m) + £(g, g) \\f\\h {m) + 2 J C(f, g)f ® g dm® 2 
for any (f,g) G T 2 and a similar expression can be given for £^ n \ 

6 Existence of a flow of maps. 

Let (S t ) t > denote the statistical solution. 

Definition 6.1 We say that (S t )t>o is a flow of maps if and only if there exists a family of 
measurable mappings (<pt)t>o from Ixfl in X such that for any f G L 2 {m) and any positive 
t, S t f = f o (f t . 
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Note that if (S t )t>o is a flow of maps, P X;W is the Dirac measure on the path {(p t (x); t > 0}. 

Definition 6.2 We say that (S t )t>o is a coalescent flow of maps if and only if (S t )t>o is 
a flow of maps and for every (x,y) G X 2 , with positive probability there exists T such that 
ip t (x) = <p t (y) for all t > T. 

let ((X t ,Y t )) t > design the two-point motion associated to the statistical solution. 

Definition 6.3 We say that (S t )t>o is diffusive without hitting if and only if (S t )t>o is not 
a flow of maps and starting from (x, x), for all positive t, X t ^ Y t . 

Definition 6.4 We say that (S t )t>o is diffusive with hitting if and only if (S t )t>o is not a 
flow of maps and (X t , Y t )t>o hits the diagonal with positive probability. 

In this section, we will give conditions under which the statistical solution is a flow of 
maps or not. 

Lemma 6.5 (S t ) t > is a flow of maps if and only if for any f G L 2 {m) and any positive t, 
E[{S t ff\ = P t p. 

Proof. It is clear that there exists Markovian kernels on X, s t (x, u, dy) such that S t f(x) = 
J f{y) s t( x ,u,dy). And s t (x,cu,dy) is the law of X t {uj') under P x ^(du)'). Asm® P-a.e, 

(S t f 2 )(x) - (S t f) 2 (x) = J (j(y)- J f(z) St (x,cu,dz)^ s t (x,cu,dy), (6.1) 

if E[(S t f) 2 ] = P t f 2 , f(f(y) - J f(z)s t (x,uj,dz)) 2 s t (x,uj,dy) = and s t (x,u,dz) is a Dirac 
measure S^x^), where ip t (x,u) is defined m <g> P-a.e. □ 

Let h G L l (m) be a positive function such that J h dm = 1. For any positive t, let 
fit be a probability on the Borel sets of X x X such that for any (/, g) G L 2 (m) x L 2 (m), 
M ® g) = f E[S t fS t g] h dm. 

Remark 6.6 (S t ) t >o is a flow of maps if and only if for all positive t, /i*(A) = 1, where 
A = {(x,x); x G X}. 

25 



Proof. If (S t )t>o is a flow of maps, there exists ip t such that S t f — f o ip t . If A and B are 
disjoints Borel sets of finite measure, 

Vk(A x B) = J E[l A (<p t (x))l B (<p t (x))] h(x) dm{x) = 0. 

This implies that /if(X x X — A) = and as /i t is a probability that /it (A) = 1. 
If /i t (A) = 1, for / E L 2 (m), fi t (f 2 ® 1 - 2f <g> / + 1 <g> / 2 ) = 0. This implies that 

/ P t / 2 hdm= [ E[(SJ) 2 ] h dm (6.2) 

and that ^'[(S'f/) 2 ] = P t f 2 - Hence (S , t ) t > is a flow of maps. □ 

Recall that we denoted by the law of the two-point motion ((X t ,Y t )) t > . 

Proposition 6.7 (S t )t>o is a flow of maps if for any positive r and any positive t, 

lim P® y) [d(X t , Y t ) >r]=0 m{dx) - a.e. 

Proof. For e > 0, let v e be the measure on X x X such that for any (/, g) e L 2 (m) x L 2 (m), 
Ve{f ® g) = J f P £ g h dm. For any (f,g) E L 2 (m) x L 2 (m), 

^A (2) (/ ® 9) = J E[SJP £ S t g] h dm. (6.3) 

As P £ S t g converges in L 2 {m <g> P) towards S t g (see remark 3.3), 

lim z/ £ P, (2) (/ ®g)= [ E[SJS t g] h dm. (6.4) 

£^0 J 

Therefore, the family of measure (^ e P/ 2 ' ) ) e >o converges weakly as e goes to towards /i t - 

3 (2) 



Assume that for any positive r and any t, lim^^^ P, , [d(X t , Y t ) > r] = 0. Let A and B 



be two disjoint Borel sets such that d(A, B) > r, then 

v £ Pt 2) (A x B) = ( f £ {x) h(x) dm(x), 
Jx 



with 



f e (x) = j P$ y) [X t E A and Y t E B] p £ (x,dy), 
where p £ (x, dy) is the kernel given by P £ . 
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As d(A, B) > r, 

fe(x) < J P$Jd(X t ,Y t ) > r] p £ (x,dy). 

For any positive f3, for m almost every x, there exists a(x) such that d(x,y) < a(x) implies 
that P$] y) [d{X t ,Y t ) >r\<(3. Note that 

fe(x) < [ p £ (x, dy) + [ P Sl)[d(X t , Y t ) > r] p £ (x, dy). 

J {d(x,y)>a(x)} J {d(x ,y)<a(x)} 

It is clear that lim e ^ f{d( x , y )>a(x)}Pz( x > d y) = ™>{dx)-a,.e. Hence, limsup / e (:r) < (3 
m(dx)-&.e and this holds for any positive (5. Therefore, lim^o fe( x ) — m(dx)-&.e and by 
dominated convergence (|/ e (x)| < 1) that 

limz/ e P t (2) (A x B) = 0. 

This implies that Ht(X x X — A) = and that (S t ) t > is a flow of maps. □ 

Proposition 6.8 // there exists a positive t, a positive r and p e]0, 1] such that for m 02 
almost every (x,y), P ( ^ > y ^{d{X t ,Y t ) > r] > p, then (S t )t>o is not a flow of maps. 

Proof. Suppose there exists a positive t, a positive r and p e]0, 1] such that for m 02 almost 
every {x,y), P™ y) [d(X t ,Y t ) > r] > p. 

Let (Bi)i(zj^ be a partition of X such that the diameter of Bi is lower than r. 

Let us suppose that /x*(A) = 1 (or that (S t ) t is a flow of maps). Then we have ^ Ht{Bi x 
Bi) = 1 and for any positive a, there exists N such that 

N 

J2MBiXBi)>l-a. 

i=i 

Since v £ Pl converges weakly towards fit, 

N 

^ ut(Bi x Bi) = 

e^O 



J^MBiXBi) = hmiy e PP(BiXBi) 

i=l 

N 

= lim V / P<>l\MX u Y t ) eB t x Bi] p £ (x,dy) h(x)dm(x) 

TV 

< lim / P$] y) [d(X t , Y t ) < r] p £ (x, dy) h(x)dm(x) <l-p 

J X x X 
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Choosing a < p, we get a contradiction. Hence /x*(A) < 1 and (S t )t>o is not a flow of maps. 
□ 

7 A one dimensional example. 

Let X — M., P t be the semigroup of the Brownian motion on M and the covariance function 
C(x,y) = sgn(x)sgn(y) (where sgn(x) denotes the sign of x with the convention sgn(O) = 1). 
Here, we have W t (x) = sgn(x)W t , where W t is a Brownian motion starting from 0. Set 
L x = sup s<t {— sgn(x)(a; + W s )} V and R x = x + W t + sgn(x)L x (it is a Brownian motion 
starting from x, reflected at 0). 

Proposition 7.1 The statistical solution S t can be written the following way 

SJ(x) = /(it?)l Lf= o + \ [f\R x t ) + /(-it?)] 1l ? >o. (7.1) 

Proof. On an extension of the probability space, it is possible to build a Brownian motion 
starting from x, X t such that W t = J* * sgn(X s )dX s (then X t is a weak solution of the S.D.E. 
dX t = sgn(X t )dW t ). Then S t f(x) = E[f \X t )\T% with T B = a(W u ;u > 0). Let us remark 
that L x t is the local time of X at and that R x = sga(x)\X t \. Set T = inf{t; L x > 0} = 
inf{t; X t = 0}. The formula (7.1) follows simply from the fact that 

E[f(X t )U> T | \X t \] = \ (f(X t ) + f(-X t )) W. □ 

8 The Lipschitz case. 

Assume X is a Riemannian manifold with injectivity radius p > 0. Let P t be the semigroup 
of a symmetric diffusion on X with generator A. Let C be a covariance inducing the metric 
(i.e with equality in (1.7)). 

We will say that C is Lipschitz if and only if there exist a positive constant k and 
< e < p such that : For any (x, y) G X 2 , with d(x, y) < e, 

A {2) d 2 {x,y) <k d 2 {x,y). (8.1) 

Remark. - d 2 (x, y) is smooth on {(x, y) G X 2 , d(x, y) < p} since p is the injectivity radius. 
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- On R d , the condition (8.1) will be checked as soon as A = \ Yli<% j<dC tj {x,x)didj + 
A k 

J2(C u (x, x) + C u (y, y) - 2C u (x, y)) < - d(x, yf (8.2) 
i=i 

and b % is a Lipschitz function for all %. 

Equation (8.2) is satisfied when C is C 2 or when C = Yla=i X a <S> X a , where X a are 
Lipschitz vector fields. In the latest case, the flow of maps can be constructed by the usual 
fixed point method for solutions of S.D.E.'s based on Gronwall's lemma. 

Let (X t , Y t ) be the two-point motion associated with the statistical solution. Set r = 
infO, d(X t ,Y t ) > e} and H t = d 2 (X tAr ,Y tAr ). 

Lemma 8.1 E™ y) (H t ) < e kt d 2 (x,y). 

Proof. By Ito's formula, 

ptAr 

H t -H = M t + / A^d 2 (X s , Y a ) ds 
Jo 

where M t is a martingale. Hence 

rtAr 

Ht — Ho < M t + kd 2 (X s ,Y s ) ds 
Jo 

< M t + k H s ds. 
Jo 

This implies that E^ y ^(H t ) — d 2 (x,y) < k J * E^ y ^(H s ) ds. Hence the lemma. □ 

Theorem 8.2 Assume (8.1) is satisfied then the statistical solution associated to P t and C 
is a flow of maps. 

Proof. Indeed, for any r < e, 

P!?Jd(X t ,Y t ) >r]< P<?Jd(X t ,Y t ) >rort>r]< ^E^(H t ) < -^d(x,y) 2 , 

which goes to as d(x,y) goes to 0. And we conclude using theorem 6.7. □ 
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9 Isotropic statistical solution on S . 

9.1 Isotropic covariance function on S d . 

On S d with d > 2, the isotropic covariance function C are given by the formula (see Raimond 



C((x,u),(y,v)) = a(t)(u,v) + P(t)(u,y)(v,x), (9.1) 
with (x,y) G S d x S d , t = (x,y) = cos<^ and (u,v) G T x S d x T 2/ 5' t2 . a and /3 are given by 

oo oo , ]_ _ +2 \ 

a(t) = ^aa^ + ^klH^-j—^t)), (9 - 2) 

i=i i=i ^ ' 

OO OO / \ 

= E^+E^f-^-^TT^)), (9-3) 

l=i i=i ^ ' 

d+1 d+1 

where ji(t) = G { \ (ty/C^ (1), Cf is a Gegenbauer polynomial, a; and bi are nonnegative 
such that Yli a i < 00 and X]; bi < oo. Using the integral form of the Gegenbauer polynomials 
(see J27| p. 496): 

TJ (cos^) = r[z(^e)] l - 1 sm d e -, (9.4) 

JO C d 

with Cd = Jq sin d 9 d9 and z((p, 6) = cos ip — i sin ip cos 6*. 

In Jl2| , it is proved that the spectrum of the Laplacian A acting on the L 2 -vector fields 
is {-1(1 + d-l), I > 1} U {-(/ + 1)(Z + tZ - 2), / > 1}. Let Q x and £>, be respectively the 
eigenspaces corresponding to the eigenvalues — 1(1 + d — 1) and — (I + l)(l + d — 2). Q t is 
constitued of gradient vector fields and X>; of divergent free vector fields. These spaces can 
be isometrically identified with the spaces Hd+i,i and Fd+i,i use d in and can be used as 



carrier spaces of the irreducible representations of 5*0(0?+ 1), T l and Q l . 

Let (ol 1 m )m and (oj 1 m )m be orthonormal basis of Qi and T>i. Then, if (z l Md )i,M and (z l MS )i t M 
are independent families of independent normalized centered Gaussian variables, 



- E E «W. +E\/d^E *UA (».5) 

i>l M 1>1 M 

is an isotropic Gaussian vector fields of covariance C given by (9.1), (9.2) and (9.3). 
Sketch of proof. The covariance of W is 

Ed Oil \ — "\ ; / ^ — "\ d bl V — "\ ; ; 

dtofi E «m * «, + £ j-^ E «» ® «- 
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Let us choose (o> 1 m )m such that a l M = ci(l, d) VE! M (where (E 1 m )m is the basis of TCd+1,1 given 
in P5] ). Then, using the fact that Hj^(p) = if M 7^ 0, for a; = and y = g 2 £> (with 
p=(0,...,0,l)), 

J2z 1 m(x)Z 1 m (v) = E ^Si^W^S^S^p) 

M M,N,K 



In [Z7p and |25[, TQ (g) is computed and it is easy from this to give the covariance of the 
gradient part of W. We can calculate the covariance of the divergent free part in a similar 
way : We choose the orthonormal basis {uj 1 m ) m of T>\ such that for M £ {!,..., d}, oj l M {p) = 



and such that for 1 < i < d, uj\{p) = 02(1, d)ei (this basis corresponds to the basis of J-d+i, 



given in [fto]). Then one have for x = gp and g G SO(d), 

d 

= J2^Mi(9)9(4(p)) = c 2 (l,d)Q l Ml (g)g(e*). (9.6) 
i=i 

Then for every (x, it) and (y, v) in TS^, 

X>y*),u><wUy),«> = (c 2 (l,d)rJ2QU9i)Q l M j (92)(9i(e i ),u)(g 2 (e^,v) (9.7) 

M 

(c 2 (/,rf)) 2 Q; i (^)(^i(e 4 )^)(^(e J )^), (9.8) 



M M 

\2r\l 



with g = g 2 gi- In [p5 |, the matrix elements Q l u(g) are calculated and it is easy from this 



to give the covariance of the divergence free part of W . □ 

Let us now introduce Sobolev spaces and related covariances. 

Let H 2 ' 3 be the Sobolev space obtained by completion of the smooth vector fields with 
respect to the norm ((—A + m 2 ) s V 1 V) 2 (with (V, V) 2 = J j | V^(a:) 1 1 2 dx), where m is positive. 
Note that the definition of H 2,s does not depend on m. 

Let a and b be nonnegative reals. Take a ; = nz^+i an d h = nzftz+i for / > 1 and 
a\ = b\ = 0. For a > 0, set G(ip) = J2i> 2 JiZT^+t1i( cos v) ■ The function G is well defined 
on [0, 7r] as |7;| < 1. 

Let Fd and F$ be real functions such that for all I > 2 

(/-l) Q+1 dim^ x F d (-l(l + d - 1)) = d (9.9) 
(/-l) Q+1 dimAxF (5 (-(/ + l)(/ + rf-2)) = d (9.10) 
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and F d (-d) = F s (-2(d - 1)) = 0. Note that when d = 2, F d = F s . 

Let II be the orthonormal projection on the space of the L 2 -gradient vector fields. 



Proposition 9.1 The covariance function defined by the sequences (a;) and (bi) is given by 
(9.1) with the functions 

a(cos<^) = aG(p) + b (cos p G(y?) + x G%)J , (9.11) 

/3(cos^) = -^-G'(<p) + b (-G{<p) + **** x G'(<p)) . (9.12) 
smyj \ (a — 1) sm<^ / 

When a and b are positive, the associated self-reproducing space is H 2 ' 2 ^ equipped with a 
different (but equivalent) norm, namely 



H = l\W\\l + \w-n)v\\l 



where \\V\\ 2 d = (F d {A)- 1 V,V) 2 and \\V\\ 2 = (F S {&)- 1 V,V) 2 . 

Proof. It is not difficult to see that the norm ||.||# given in the proposition is the norm on 
the self-reproducing space associated to C . 
Now since (see [0) 

dim ^ = (d-t)!(7-l)! (<f + 2f " 3)(d+1) ' 
, m (d + l-3)\ , ,d(d+l) 

for A — > oo, A^F rf (A) = 0(1) and X £ ^Fs(X) = 0(1). This implies that ||.||# and the 
norm used to define H 2 '~5~ are equivalent (when a and b are positive). And we get that the 
self-reproducing space associated to C is H 2 ' 2 ^ . □ 

Remark 9.2 If a or b vanishes, the self-reproducing space is H 2,9 ¥ restricted to divergent 
free vector fields or gradient vector fields. 

9.2 Phase transitions for the Sobolev statistical solution. 

Let P t be the semigroup of the Brownian motion of variance (a + b)G(0) and St be the 
statistical solution associated to P t and C. 
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Let (X t ,Y t ) be the two-point motion. Let ipt = d(X t ,Y t ). Since h(x,y) = d 2 (x,y) is a 
C 2 -function, h belongs to V(A^) and since X t and Y t are solutions of an S.D.E. like (3.3), 

is a diffusion on [0, n 2 \ and is solution of an S.D.E. ip t is also a diffusion on [0, n] (note 
that d(x,y) a priori does not belong to T)(A^)). This diffusion is eventually reflected (or 
absorbed) in and it. Its generator is L = o~ 2 ((f)-^2 + Kv 9 )^ ( see R-aimond f25|), with 

cr 2 (<y?) = ot(l) — a(cos(p) cos (f + [3(cos(p) sin 2 cp, (9.13) 

Utp) = ^ — — (a( 1) cos (p - a( cos ip)). (9.14) 
sm if 

The generator of ip 2 is V = a 2 {x)-^ + K x )^^ with 

a 2 (x) = 4xa 2 (Vx), (9.15) 
b(x) = 2a 2 (v / x) +2v^6(v^). (9.16) 



Lemma 9.3 If a > 2, the statistical solution is a flow of maps. 

Proof. We have A^d 2 (x, y) = 2a 2 (d(x, y)) + 2b(d(x, y))d{x, y). When a > 2, then G is C 2 , 
this implies that a is C 2 and /? is continuous. Hence equation (8.1) can be checked. □ 

Suppose a + b > and let r\ = 
Theorem 9.4 For any a e]0, 2[, 

• For d = 2 or 3 and rj < 1 — \, the statistical solution is a coalescent flow of maps. 

• For d = 2 or 3 and 1 — ^<r/<| — ^~ 2 ^ , £/ie statistical solution is diffusive with 
hitting. 

• For d = 2 or 3 and rj > | — ( d ~ 2 ) or /or d > 4, i/ie statistical solution is diffusive 
without hitting. 

Remark. The same phase transition appears in the M. d case (see theorem 10.1 below). It 
has been independently observed, in the context of the advection of a passive scalar, by 



Gawedzky and Vergassola 14 
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Lemma 9.5 For a e]0,2[, we have 

• G is differentiable on ]0, ir[. 

• lim^ 0+ = i; f? ^ht^ sin^ ^ = KG(0). 

• lim^ 0+ = -a /; pg^ 1 sin<* (9 j*^ = -a* G(0). 

The proof of lemma 9.5 is in appendix A. From this lemma, we get as ip goes to 

«(cos ip) = (a + b)G(0) - (a + (l + -^-^J b^j KG{0)ip a + o{ip a ) (9.17) 
/3(cosip) = a (a- — — j KG(0)ip a ~ 2 + o(ip a ~ 2 ). (9.18) 



d- 1 



Hence, 



ct 2 (<^) = ( a + 6)irG(0)(a + l-ar/)^ Q (l + o(l)) (9.19) 
b{ip) = (a + 6)KG(0)(rf-l + ar7)^ a - 1 (l + o(l)). (9.20) 

In order to prove theorem 9.4, we need to study the two-point motion. Because of 
isotropy, it is enough to study the diffusion tp t . This diffusion satisfies an S.D.E. until it 
exits ]0, 7r[. 

Let s be the scale function of the diffusion ip t , 

'» b(ip) 



s(x) = / exp 



X 



dip 



dy, with (xq,x) g]0, ir[ . 



/, ^) 

Let x G {0, 7r} and T x = inf {£ > 0; ^ = x}. Using Breiman's terminology (see [Q p. 368- 
369), x is an open boundary point if T x = oo and is a closed boundary point if T x < oo. 
Note that x is an open boundary point if \s(x)\ = oo . 

Firstly we are going to show that 7r is an open boundary point. Then : 

• When d = 2 or 3 and rj < 1 — -4r, we prove that is an exit boundary point (this 
implies that the statistical solution is a coalescent flow of maps). 

• When d = 2 or 3 and 1 — ^|-<?7<| — ^ d ~^ , we prove that is an instantaneously 
reflecting regular boundary point (this implies that the statistical solution is diffusive 
with hitting). 
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When rj > ~ — ^ 2 ^ , we prove that is an open entrance boundary point (this implies 
that the statistical solution is diffusive without hitting). 



Lemma 9.6 n is an open boundary point. 

Proof. It is easy to check that s(tt— ) = oo using the fact that a(l) + a(— 1) > : 

a(l) + a(-l) = (a + b)G{0) + (a - 6)G(tt) > (a + 6)G(tt) + (a - 6)G(tt) > 0. □ 

Since 7r is an open boundary point, we now study the behavior of ip t at and near 0. 

Lemma 9.7 Ifr]>\ - s(0+) = -oo and if 7] < \ - s(0+) > -oo. 

Proof. Let us note /i = ^~ . Then, we have that Jr^y = + an d f° r an y positive 
£ there exist positive constants C\ and C2 such that for y < xq, 

r m 
L ^) 



Ciy-» +£ < exp 



dip 



<C 2 y 



—/l—E 



(9.21) 



From this, we see that s(0+) = —00 if \i > 1 (or if 77 > | 



1 _ (rf-2) ■ 



2a 



and s(0+) is finite if \i < 1 



(or if 77 < 



1 (d-2) 



2a 



□ 



Lemma 9.6 and 9.7 implies that (see theorem VI-3.1 in |16[) if 77 < 



1 (d-2) 



we have 



T < 00, T n = 00 a.s. and if 77 > ~ — ^ 2 ' , is an open boundary point and we have 
liminf ip t = and limsup^ = ir a.s ("0* is recurrent). 

Remark 9.8 When d > 4 and a g]0, 2[ ; ~ — *- d ~ 2 ^ < 0. 27ms implies that liminf ip t = and 
limsup ipt = n a.s. 

Since n is an open boundary point, ip t £ [0, 7r[ for every positive t and ^ is a solution of 
the S.D.E. 

dip] = V2a(^) dB t + b(ip t ) dt. (9.22) 

Note that is a solution of this S.D.E. (since a(0) = 6(0) = 0). The solutions of this S.D.E. 
might be not unique. 

Let m(dx) be the speed measure of the diffusion : 



m(dx) = l]o, n [(x) exp 
with xq G]0, 7r[. 



^feO 



dp 



dx 



+ m({0})5 = c/(x) da; + m({0})5 , 
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Lemma 9.9 If-q>\- is an entrance open boundary point. 

Fro of. When 77 > f - ^1, is an open boundary point. From Proposition 16.45 in M, 
is an entrance boundary point if and only if J 0+ \s(x)\m(dx) < oo. For any positive e, there 
exists a positive constant D such that, for any x g]0, x [, 

\s(x)g(x)\ < D x (f-°-)AO x -f- £ +i < D x i-«-2 £ _ 

This shows that J Q s(x)m(dx) < oo (choose e such that 2e < 2 — a). □ 

This lemma implies that when rj > ~ — , there exists a positive t, a positive a and 
p G]0, 1[ such that for any x g]0, tt[, P x [ip t > a] > p. Proposition 6.8 implies that S t is not a 
flow of maps and since is open, St is diffusive without hitting. 

Let now d G {2, 3} (when d > 4 we always have rj > | — ) . 

Lemma 9.10 If rj < ~ — > is a closed boundary point. 

Proof. From Proposition 16.43 p. 366 in M, T is finite or the boundary point is closed if 
and only if for any b g]0, 7r[, f^ \s(x) — s(0)\m(dx) is finite. 

We have \s(x) — s(0)\g(x) ~ Jq^P 
\s(x) — s(0)\g(x) = 0(x 1 ~ a ). This implies that J " \s(x) — s(0)\m(dx) is finite. This proves 
that T is finite a.s. □ 

Lemma 9.11 If -q < 1 - 4, is an exit boundary point. 

Proof. In |D, is an exit boundary point if and only if m(]0,x[) = oo for all x G]0,7r[. 
This is the case if fi — a < — 1 (or if r] < 1 — %). Note that for d = 2 or 3 and a g]0, 2[, 

i _ 4 < i _ D 

Lemma 9.11 implies that when r] < 1 — the diffusion ^ is absorbed at 0, and for any 
positive r, 

hm Pg> y) [d(X t ,Y t ) >r} = \imP^ t > r] = 0. 

Now, applying proposition 6.7, we prove that the statistical solution is a flow of maps and 
this is a coalescent flow of maps (since is an exit boundary point). 
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L 



dip 



exp 



P 

Jx, 



dip dy. Hence 



, is a regular boundary point. 

Proof. In HJ], we see that is regular if m(]0,x[) < oo for all x G]0,7r[, which is the case 

. □ 

, the two-point motion hits the diagonal. But there is no 
uniqueness of the solution of the S.D.E. satisfied by ip t since might be absorbing or (slowly 
or instantaneously) reflecting. In order to finish the proof of theorem 9.4, we are going to 
prove that is instantaneously reflecting. 

To prove this, for e G]0, 1[, let us introduce the covariance C £ = (1 — e) 2 C (then, if W t 
is the cylindrical Brownian motion associated to C, (1 — e)Wt is the cylindrical Brownian 
motion associated to C e ) and Sf be the statistical solution associated to P% and C £ . 

For / G L 2 (dx), Sff = J2 n>0 J™' 6 f, where J™' £ / is the nth chaos in the chaos expansion 
of Sff 0. It is easy to see that J™' 6 f = (1 — e) n J™f, where J™f is the nth chaos in the chaos 
expansion of Stf, hence 

E[(Stf - S t f) 2 } = £(1 - (1 - e) 2n )E[(J?f) 2 }. (9-23) 

n>l 

Hence it is clear that the L 2 (P)-limit as e goes to of Sff is S t f. 

Let (X £ t ,Y t £ ) be the Markov process associated to P t (2) ' e = E[S £ t m ] and tp £ = d(X £ ,Y £ ). 
ijj £ is a diffusion with generator L e . It is easy to see that L e = (1 — (1— e) 2 )Li + (l — e) 2 L (note 
that Ai 2) = A®I+I®A+(l-s) 2 C = 4 2) + (l-e) 2 (A( 2 )-4 2) )), andL £ = a 2 ^)£,+b £ (ip)^, 
with 

a 2 M = (l-(l- £ )> 2 H + (l-e)VH, (9.24) 
b £ (<p) = (l-(l- £ ) 2 )M^)|(l- £ ) 2 %). (9.25) 

Let us remark that L\ is the generator of the diffusion distance between two independent 
Brownian motions on S d . Note that as (p goes to 0, 

a 2 A V ) ~ a 2 (0) = 2(a + b)KG(0) and b x Up) ~ ^—^-(a + b)KG(0) (9.26) 
<P 

2 Note that Sf = Qi oe .(i- s )St, where Q a is the Ornstein-Uhlenbeck operator on the Wiener space (used in 
Malliavin calculus : see p4[). 



Lemma 9.12 If r] G 



_d_ 1 

a 2 ' 2 



(d-2) 
2a 



when r] G 



1 - 



d 



i.d-2) 
2a 



When 7] G 



i _ d, 1 

J- : 1 r, 



(d-2) 
2a 
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wdo%(<p) = (l-(l-e) 2 )a2(v9)(l+0(v9 Q ))and6 £ (v9) = (l-(l-e) 2 )6i( ¥ »)(l+0( ¥ » a ). Studying 
the scale function s e of ip £ , we get that s e (0+) = Si(0+) = — oo (as two independent Brownian 
motions cannot meet each other on S d ). We still have s £ (ti— ) = oo. Hence ipf e]0, 7r[ for all 
positive t. 

Let m £ be the speed measure of Let g £ (x) = m £ (dx)/dx. As m e (]0, 7r[) < oo, m £ is 
an invariant finite measure for the diffusion ip £ . As lim e ^ = °" 2 and lim e _>o b £ = b, we get 
that lim £ ^ g £ (x) = g(x). Let us note e' — 1 — (1 — e) 2 and let 

,(£ '^P*)TF7KM' (9 - 27) 

This function increases with e' if > As — ~ (rf — 1 — a)- as <z> goes to 

and as (d — 1 — //) is positive, there exists ip such that for any (p < </? , /(e', y?) > ^(^y 
/(0,^>) and for e' < 1/2, 



and as (d — 1 — ji) is positive, there exists ip such that for any (p < </? , /(e', y?) > T^fl = 



2 ( /*v^o ^ 



where C vo = sup ee [ 0jl ] exp ( /(e', </?) dy>) < oo. The Lebesgue's dominated convergence 
theorem implies that g £ converges in -^ 1 ([0, n]) towards g. 

Let / and g be continuous functions, then E[f(X £ )g(Y t £ )} = E[Sff(x)Sfg(y)]. Since Sff 
and S £ g converge respectively towards S t f and S t g when e goes to in L 2 (P), we get that 
(X s , Y t £ ) converges in distribution towards (X t , Y t ) when e goes to 0. This also implies that 
tp £ converges in distribution towards ip t when e goes to 0. 

Since m e is an invariant measure, for any continuous function / on [0, n], we have 

J E[f(r t M £ = x] m £ (dx) = J f dm £ . (9.28) 

Since 

j E[f(r t Wo = A m £ (dx) - J E[f(ij t )\ij = x] g(x)dx < 

< ll/IU f \9e(x) - g(x)\ dx + r (E[f(r t M = x]- E[f(^ = x}) g(x)dx , 
Jo Jo 

we get that (because g £ converges in L 1 ([0,7r]) towards g and ip £ converges in distribution 
toward tp t -) 

J E[f (ip t )\i) = x] m(dx) = lim J E[f(ijj £ )\ipl = x] m £ (dx) = 1 ™J f dm e = J f dm - 
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This implies that g(x)dx is an invariant measure for %p t and m(dx) = g(x)dx. Since 
m(]0, x[) < oo for all x e]0, ir[, the diffusion ip t is not absorbed in and is reflected in 
0. 

In this case, is a closed regular boundary point. This point is instantaneously reflecting 
since m({0}) = 0. This implies the existence of a positive t, a positive r and p g]0, 1] such 
that for any x g]0, 7r[, P x [^t > r ] > V- Then, applying proposition 6.8, the statistical solution 
is not a flow of maps. This finishes the proof of theorem 9.4. □ 

For a > 2, the statistical solution is an isotropic Brownian flow of diffeomorphisms. In 
Raimond p5| , the Lyapunov exponents of this flow are computed. The sign of the first 
Lyapunov exponent X\(a,d) describes the stability of the flow. It is unstable if Ai > and 
stable if Ai < 0. The computation of \i(a,d) gives 

(d-A)a + db^., fd-l\ tfl „ 1>/ . ,/2(d - I) a + db\ . . 

Al = d + 2 C(Q " 1} + (—2) [id ~ A)a + db]Cia) - d { dU ) C(a + 1} ' 

(9.29) 

where C( a ) = Yli>i h ls the zeta function. Therefore, we have Ai(a, d) = if and only if 

-{d - 4)C(a - 1) - (d - l)(d - 4)C(a) + 2rf(rf - l)C(a + 1) 
^ ^ ' ' 4C(a-l) + 4(d-l)C(a) + d(d-2)C(a + l) ' 1 J 

It is easy to see that for fixed 77, lim Q ^ 2 + Ai(a, d) = +00 if d > 4 or if 77 > ~ — = -^p and 
that lim Q ^ 2 + Ai(a, d) = — 00 if 77 < Remark that lim Q ^ 2 - 1 — ^ = lim«^2- § — 2 = 
This shows that coalescence appears when Ai goes to —00 and splitting appears when 
Ai goes to +00. 

The results of this section is given by phase diagrams in appendix B. 

10 Isotropic statistical solution on M, d . 

10.1 Stationary and isotropic covariance functions on M. d . 



On M. d with d > 2, the stationary isotropic covariance function C are (see Le Jan |20|) such 
that C ij {x, y) = C ij {x - y), for (x, y) G R d x R d , with 

C«(z) = S*B N {\\z\\) + ^(B L (\\z\\) - B N (\\z\\)), (10.1) 

\\ z \\ 
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with 

BAr 



J j cos(pu 1 r)u 1 uj(du)(F L (dp) - F N {dp)) + J j cos(pu 1 r)u(du)F N (dp)(10.2) 

Bn{ t ) = J J cos(pu 1 r)u 2 ,uj(du)(F L (dp) — F N (dp)) + J J cos(puir)ui(du)F N (dp)(10.3) 

F L and F N being finite positive measures on IR + . uj{du) is the normalized Lebesgue measure 
on S d ~ l . Fl and Fn represent respectively the gradient part and the zero divergence part 
of the associated Gaussian vector field. 



„d-l 



For a and m positive reals, let F(dp) = e — j+^dp, F L (dp) = aF(dp) and F N (dp) = 

(p 2 +m 2 ) 2 

-^■F(dp), where a and 6 are nonnegative. In the Fourier representation, (c is a positive 
constant) 

CHk) - 4W + (-H + ^ (% - g)) ■ dO.4) 

Notice that in the Fourier representation, the Laplace operator on vector fields is given by 
the multiplication by — ||A;|| 2 and the projection n on gradient vector fields (in the L 2 space) 
by 0. (i.e if V is a vector field and V\k) its Fourier transform, (irV)\k) = ^V j (k).) 

Therefore, given a L 2 vector field, U^(y) = j E« C^(x — y)V l (x) dx can be expressed as 
c(-A + m 2 Y^[aW + ^(J - ir)V). Since ([/, U) H = (U, V) 2 = j(U{x), V{x)) dx, the 
self-reproducing space appears to be the L 2 -Sobolev space of order s = (defined the 
same way as in secion 9.1) equipped with the norm 

ll^|| 2 = i||^|| 2 + ^||(/- 7 r)y|| 2 , 

where 

\\V\\ 2 s = - c ((-A + m 2 yv,V) 2 . 

Note that if a or b vanishes, the self-reproducing space is H 2,9 ¥ restricted to divergence 
free vector fields or gradient vector fields. 

10.2 Phase transitions for the Sobolev statistical solution. 

Let P t be the semigroup of a Brownian motion on R d with variance (a + 6)F(M + ). Let S t 
be the statistical solution associated to P t and C. If a > 2, C is C 2 . Hence equation (8.2) 
is satisfied and the statistical solution S t is a flow of maps. 
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Suppose a + b > and let r\ = Then we have the theorem. 

Theorem 10.1 For any a e]0, 2[, 

• For d = 2 or 3 and r] < 1 — , the statistical solution is a coalescent flow of maps. 

• For d = 2 or 3 and 1 — ^<n<| — ^~ 2 '* ; the statistical solution is diffusive with 
hitting. 

• For d = 2 or 3 and i] > | — ^ d ~ 2 ^ or /or d > 4, £ne statistical solution is diffusive 
without hitting. 

Remark. The results of this theorem are exactly the same as for the sphere. 

Proof. Let us study the two-point motion (X t ,Y t ) starting from (x,y) (with x ^ y). 
Then r t = d(X t ,Y t ) is a diffusion in M + (eventually reflected in 0), with generator L = 
+ b{r)f r (see Le Jan @), with 



a\r) = B-B L (r), (10.5) 
b(r) = (d-l f- BN{r \ (10.6) 



where B = B L {0) = B N {0) = s±*F(R+). 

Lemma 10.2 For a g]0,2[, as r goes to 0, 

i) J J cos{pu 1 r)u(du)F(dp) = F(R + ) - a x r a + o(r a ). 

ii) J J cos{puir)u\uj(du)F(dp) = F ^ ^ — a 2 r a + o(r a ). 
Hi) J J cos(puir)u\oj{du)F (dp) = — a 3 r a + o(r a ). 



with ct2 = TT^ai, a 3 = -77 — tti and 



Oil = c d 



l-cosxl-^rW [ * (cos 9) a (sin 9) d - 2 d9 



o x / Wo 
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Proof. For r > 0, making the change of variable x = puir, 

p p pi pco 

/ / (1 — cos(pu 1 r))u(du)F(dp) — c d / / (1 — cos(pu 1 r))(l — u^) 



d-2 , p d 1 dp 
2 dux 



(p 2 + m 2 ) a" 



= r a c d 



1 / />oo 



o \Jo 



(l-coss) - X f,.^ ) «?(l-u?) d » 2 d Ul . 
(ar + r l u\m l ) 2 



As lim r ^ / °°(1 — cos a;) 



(x 2 +r 2 ufm 2 ) 



-^dx = / °°(1 - cos:t)^#ft < oo, we get that 



a) = ai, 



with J(n,t) = / o 2 (cos0) t (siii0) n d6 = ±B(2±i,*±I) for t > and n G N, and = 
F r(x+y) • ^his shows i). ii) and iii) can be obtained the same way with 

f°° dx 
oti — c d / (1 - cosa;) — — rl(d -2, a + 2) 



x 



and ai = a 2 + (d — l)a 3 (note that J u\ui{du) — |). It is easy to see that for ct > and 
d> 1, 

ct + 1 

J(d - 2, a + 2) = -r^I(d -2, a). 

di ~\~ cx 

Therefore, 0:2 = With the relation ai = a 2 + (d — 1)0:3, we get that 03 = ^^ai- U 

Remark 10.3 As z goes to ; 



C^'(z) = B5 ij - 



Oil 
d-1 



((d -l)a + (d + a- l)b)5 ij - a{{d - l)a - b) 



1?3 



z z- 



\z\\ a (l + o(l)), 



Let us note that the dependence on m only appears in B . 



From this lemma, it is easy to see that as r goes to 0, 

a 2 (r) = ( ° + 6)ai (a + l-aiy)r a (l + o(l)), 



6(r) = 



d + a 
(a + b)a>\ 



{d-l + ar])r a -\l + o{l)). 



d + a 

Note that we get the same behaviour of a and b around as in section 9.2. 
As in section 9.2, let us study s, the scale function of the diffusion r t . 



(10.7) 
(10.8) 
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Since Bl(t) and Bn(t) converge towards as r goes to oo (as Fourier transforms of finite 
measures), we get that as r goes to oo, log(s'(r)) ~ (1 — d) log(r). Therefore s(+oo) is finite 
if and only if d > 3. 

We also see that s(0+) = — oo if rj > | — ^ d ~ 2 ' and s(0+) is finite if rj < ~ — ^ d ~^ . 

Let m be the speed measure of the diffusion. Let us study the boundary point 0. 

As m(]0, x[) < — oo for any positive x if 77 > 1 — as in section 9.2, with a similar proof, 
we can prove that if -q e]l - | - the diffusion r t is instantaneously reflecting at 0. 

The only thing there is to change in the proof is to take the test function / in (9.28) with 
compact support and to remark that g £ converges towards g in L\ 0C (JH + ). 

If 77 < 1 — 4r (note that 1 — -4r < § — ^-^), is an exit boundary point and the diffusion 
is absorbed by 0. 

Therefore, we get that 

• If d > 3 and 77 e]1 — ^, | — ^ [, T t is instantaneously reflecting at and is transient. 
In this case, as in section 9.2, (S t )t>o is diffusive with hitting. 

• If d = 2 and 77 e]l — ^, | — [, r 4 is instantaneously reflecting at and is recurrent. 
In this case, as in section 9.2, (5 , t ) t > is diffusive with hitting. 

• If d > 3 and rj < 1 — r< is absorbed at with probability an d converges 
towards +00 with probability ^[^-^(oj • ■'■ n ^ ms case ! as m the section 9.2, (St)t>o is a 
coalescent flow of maps. 

• If d = 2 and 77 < 1 — ^, r t is absorbed at a.s. In this case, as in section 9.2, (S t )t>o 
is a coalescent flow of maps. 

If 77 > I — then we have that s(0) = —00. In this case, is an entrance boundary 
point as J 0+ \s(x)\dm(x) < 00. r t is recurrent if d = 2 and transient if d > 3. As in the 
section 9.2, we prove that (S t )t>o is diffusive without hitting. □ 

For a > 2, the statistical solution is a stationary isotropic Brownian flow of diffeomor- 
phisms. In Le Jan PDA, the Lyapunov exponents of this flow are computed. The sign of the 
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first Lyapunov exponent Ai(a, d) describes the stability of the flow. It is unstable if Ai > 
and stable if Ai < 0. The computation of Xi(a,d) gives (see EC 



Ai= ((d-4)a + db) I p 2 F(dp), (10.9) 



2(d + 2 

Therefore, we have Ai(a, d) = if and only if d < 4 and 

ri = V (d) = (10.10) 

As in the section 9.2, we see that for fixed 7), lim Q ^2+ Ai(a, d) = +oo if d > 4 or if 7] > 
\ — = and that lim Q ,_ > 2+ Ai (a, d) = — oo if r\ < -M. This shows that coalescence 
appears when Ai goes to — oo and splitting appears when Ai goes to +oo. 
Remark that lim a _ 2 _ 1 - 4^ = lim a ^ 2 - \ - = 
The results of this section are given by phase diagrams in appendix B. 

11 Reflecting flows. 

Let D be an open convex domain in M. d with C 1 boundary dD. Let d be the Euclidean 
metric in M. d . For any x G dD, we denote n(x) the directed inward unit normal vector to 
dD. 

Let P t be the semigroup of the Brownian motion in D reflected on dD. P t is associated 
to the Dirichlet form (£,.F), where T = H\D) = {/ e L 2 (D,dx), |V/| G L 2 (D,dx)} 
equipped with the form | |V/| 2 dx (see [p]] , 1.3.2). Let C(x,y) be a covariance function 



in D x Z) such that C l i(x,x) = 8^ and satisfying (8.1). 

We can construct a statistical solution associated to Pt and C. Let P t be the semigroup 
of the two-point motion (JQ, 5^). Let P^ 2 \ be the law of the two-point motion. 

We know that X t and Y t are two diffusions in D reflected on dD. Let (ft and i/j t denote 
the local times of X t and Y t on dD. 

Lemma 11.1 For h(x,y) = d 2 (x,y), P^h(x,y) < h(x,y)e ct . 
Proof. Let us note 



= A x + A y + J2c iJ (x,y)d Xt d x 



i,0 
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From (8.1) and the Lipschitz conditions, we get that 

L^h(x,y)<C h(x,y). 

Using Tanaka's formula, there exists a martingale M t such that 

h(X t ,Y t )-h(x,y) = M t + f L^h(X s ,Y s )ds (11.1) 

Jo 

+ / (V x h(Xs,Y s ),n(X s ))d<p s + [ {V v h(X s ,Y a ),n(Ys))d^ s {11.2) 
Jo Jo 

As V x h(x, y) = 2(x — y), using the fact that D is convex, we get that for x E dD 

(V x h(x,y),n(x)) < 0. 

This implies that 

h(X t , Y t ) - h(x, y)<M t + C [ h(X s , Y s ) ds. 

Jo 

Taking the expectation, we get that P^h(x,y) — h(x,y) < C J* h(x,y) ds. Hence the 
lemma. □ 

Theorem 11.2 The statistical solution is a flow of maps. 
Proof. This is the same proof as the proof of theorem 8.2. □ 

A Proof of lemma 9.5. 



Take G]0,7r[. At first, we are going to prove that /(</?) = J2i>2 (PTp+i i^7'( cos V 9 ) 
finite. A S i = /~e-'» S <»- 1 fe, 

Jo Jo frf \z{(f,0)\ L{a)c d 



i>i 

TV pQO 



< / / UM dsd6 = 2 / f Vid (s) ds d6, (A.2) 
Jo Jo Jo Jo 



with f^e(s) = tJ-^fl), r s (a)cd ■ It is easy to see that 



1 f 00 e- s s a ~ l 



U,e( s ) ds < tv \ / 71 z^ ds < oo. (A.3) 

i T(a)c d J 1 (1-e s ) 
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On the other hand, 



f M s ) ds ^ wi^T I i ^H7Z mi = wfvr^W- ( A - 4 ) 



1 Z" 1 rfs 1 

r(a)c d> / 1 - e- s |z(</?,0)| r(a)c d ' 

Let x v (0) = - log \z(ip, 0)\, then F v {9) = f^$ +1 -^- t . As lim^ + = °> we havc 
Fip(9) ~ — log 2^(0) as 6 1 goes to 0. From this, we see that F v (9) = 0(\og9) as 9 goes to 0. 
This implies that I(<p) is finite. 

Now, applying the derivation under the integral theorem, we prove that G is different iable 
on ]0, 7r[ and that for ip g]0, 7r[, 

G'(<f) = J2L la sin9 T d (A - 5) 

Er -s s a \^~k Z ( {p,e ^ a-i ■ dn dsd() (\ a\ 

[e Z ^ 6)]J£ zJ^V S Smd f(a^ (A - 6) 



1>1 "0 
<*7r /"OO 




JO 



1>1 



/o Jo l-e-°z(f,9) S Sm9 T(a)c d ^ 



As z(p, it — 9) = z(ip, 9), 



Jo Jo b(s,p) + cos 2 9 snap T(a)c d 
with a(s, ip) = e _4osy and 6 ( s ' V 9 ) = (1 ~-2 3 3 s -°2^ )2 . Changing of variables (s = tip), 
G'(p) _ [°° a(tp, p) - sin 2 9 ^ pcosp ^, ^ fl 2rft rffl 



</? Q 1 Jo Jo b(tip, ip) + cos 2 9 sin 99 r(a)cd 

I(t,<p,9) dt d9. (A.9) 




'0 jo 

Let e > 0, there exists a positive constant C £ such that for any t G [0, e], 

< I(t,p,9) < C £ t a -\ (A.10) 

Remark also that 

I(t,f,9)<C d , a ^f^ )2 xt"-\ (A.ll) 

2 —x 

where C d)Ct is a positive constant. Let C = C d)Cc sup ;E>0 ^_^- x y < oo, then, for any positive t 

0</(t,^^)<Cr 3 . (A.12) 
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As F(t) = Cet^loKtKe+C t a - 3 l t>£ belongs to ^(dO^dt) for a e]0, 2[, lim^ a(^, V 5 ) = 
1 and lim^o </?) = by the Lebesgue dominated convergence theorem, 

G'((p) H f°° cos 2 9 a , . dn 2d6dt T ^ 
lim — = - / / — * a-1 sui d — . = -aAT. A.13 

We have proved the second limit. The first limit is easy to obtain as 

G(Q)-G(<p) = - f G\x) dx 
Jo 

= -K V a + o{^ a ). 

This finishes the proof of the lemma. □ 
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Phase diagrams for the Sobolev statistical solutions. 
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i i i i 


diffusive without hitting 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
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Figure 1 : Phase diagram on S 2 . 


Plane 6=2 


stable flow 


unstable flow 






diffusive with hitting 


diffusive without hitting 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
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Figure 3 : Phase diagram on S* 3 . 
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Figure 6 : Phase diagram on S* 5 . 
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-—-""^ flow 
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diffusive without hitting 
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1 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 7 : Phase diagram on S 50 . 

Let us remark that when a < 2, the diagrams are exactely the same for the sphere and for 
the plane. For the sphere, we see that, for a > 2 and r\ < 2 — ' ^ ne ^ ow S e ^s stable 
when d goes to oo : (9.30) implies that \im d ^ OQ r](a,d) = 2 — for a > 2. We see 

that, for any d and 77 G [0, 1[, the flow gets stable when a goes to 00 : (9.30) implies that 
lim^oo 7](a, d) = I. 
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